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ABSTRACT: Strangeness production was investigated during the equilibration of 
a gluon dominated parton plasma produced at RHIC and LHC energies. The time 
evolution of parton densities are followed by a set of rate equations in a 1-dimensional 
expanding system. The strangeness production will depend on the initial chemical 
equilibration level and in our case the parton densities will remain far from the full 
equilibrium. We investigate the influence of gluon fragmentation on final strangeness 
content. 



Introduction 

Strangeness production was widely investigated during last decade as one of 
the possible quark-gluon plasma signatures . To refresh our knowledge about 
the time dependence of strangeness content in case of parton gas, let us con- 
sider ultrarelativistic heavy ion collisions in the framework of pQCD-inspired 
parton interaction models |2|, |3|, f|]. In this framework, hard or semihard scat- 
terings among partons dominate the reaction dynamics. They can liberate 
partons from the individual nucleons, thus producing large amount of trans- 
verse energy in the central region ||, and drive the initially produced parton 
system toward equilibrium [||, |, [?[]. 

Roughly speaking, ultrarelativistic heavy-ion collisions in a partonic pic- 
ture can be divided into three stages: (1) During the early stage, hard or 
semihard parton scatterings, which happen on a time scale of about 0.2 fm/c, 
produce a hot and undersaturated parton gas. This parton gas is dominated 
by gluons and its quark content is far below the chemical equilibrium value. 
Interference and parton fusion also leads to the depletion of small x partons in 
the effective parton distributions inside a nucleus, reducing the initial parton 

1 The talk was presented by P. Levai on S'95 Conference, Tucson, Jan. 3-7 1995. E- 
mail : plevai@sunserv.kfki . hu 
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production ||. (2) After two beams of partons pass through each other, the 
produced parton gas in the central rapidity region starts its evolution toward 
(kinetic) thermalization mainly through elastic scatterings and expansion. The 
kinematic separation of partons through free-streaming gives an estimate of 
the time scale r iso ~ 0.5 — 0.7 fm/c ]7| [|, when local isotropy in momentum 
distributions is reached. (3) Further evolution of the parton gas toward a fully 
(chemically) equilibrated parton plasma is dictated by the parton proliferation 
through induced radiation and gluon fusion. This chemical equilibration can 
be described by a set of rate equations for gluons and quarks with different 
flavours. Due to the consumption of energy by the additional parton produc- 
tion, the effective temperature of the parton plasma cools down considerably 
faster than the ideal Bjorken's scaling solution. Therefore, the life time of 
the plasma is reduced before the temperature drops below the QCD phase 
transition temperature [^]. 

Strangeness production could also be divided into the above three different 
contributions in the history of the evolution of the parton system. However, 
here we will skip the very complicated description of steps (1) and (2) and 
simply we consider the result of the HIJING Monte Carlo code for these early 
stages. We will concentrate on strangeness production in stage (3), investigat- 
ing the strangeness evolution and saturation. We describe the equilibration of 
the initially produced hot and undersaturated parton gas following Ref. 0. 
For gluons we apply the result of an improved perturbative QCD analysis of 
Landau- Pomeranchuk-Migdal effect (see Ref. and Ref. fTTH ). 



In our analysis we include one more step of strangeness production, namely 
the gluon fragmentation. We assume a two step hadronization process at 
a certain critical temperature: (a) gluons decay into quark-antiquark pairs 
forming quark-matter from the parton gas; (b) the quarks and antiquarks will 
be redistributed forming colorless hadrons. For simplicity, we will use the 
well-known string breaking factors for estimating the different flavour yields 
in the gluon decay reaction (g — > qi + gj. However more work will be 
needed to find out the medium- dependent values of these factors. The hadron 
formation by quark-antiquark redistribution is beyond of the scope of our 
recent investigation. However, the strange/non-strange quark-pair ratio will 
not change in this type of hadronization. One can find a detailed description 
of such a hadronization mechanism in Ref. ||12|| . 

We will study the strangeness production during the parton evolution. We 
also discuss the influence of the uncertainties in the initial parton density 
on the final saturation level of different flavours. We investigate the extra 
strangeness yield of the gluon fragmentation and we compare the change of 
the strange/non-strange quark ratio in a chemically non-equilibrated parton 
system to a fully equilibrated one. 
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Initial conditions: hot and undersaturated 
parton gas 

The initial conditions for the time evolution of thermalized but chemically 
non-equilibrated parton gas can be obtained from one of the microscopical 
models. Currently there are many models for incorporating hard and semi- 
hard processes in hadronic and nuclear collisions 0, |[]. Here we will use 
the result of the HIJING Monte Carlo model || to estimate the initial parton 
production at time, r iso . 

Since we are primarily interested in the chemical equilibration of the parton 
gas which has already reached local isotropy in momentum space, we shall 
assume that a parton distribution can be approximated by thermal phase 
space distribution with non-equilibrium fugacity A: 

f(k;T,\) = \(e u - k/T ±\Y\ (1) 



where w M is the four-velocity of the local comoving reference frame. When 
a parton fugacity A is much less than unity as may happen during the early 
evolution of the parton system, we can neglect the quantum corrections in 
Eq. ([I]) and write the momentum distributions in the factorized form, 

f(k;T,X) = X(e^±l)~\ (2) 

Using this distribution in the comoving frame for a multi-component parton 
gas, one can obtain the total parton density, energy density and pressure as 

n = n g + J2(n qi +ngJ = (\ g a 1 + Y J ^i b i( x i))- T3 ( 3 ) 

i i 

£ = ^ + E(^ + ^) = (V2 + E A ^(x i ))-T 4 (4) 

i i 

P = P g + J2( P i* + P gJ = 77(-V3 + E kb2(xi)b 3 (xi)) ■ T 4 . (5) 
i 6 i 

Here a\ = 16C(3)/7r 2 « 1.95 and a 2 = a 3 = 87r 2 /15 « 5.26 for the Bose 
distribution of gluons. For fermions (i = u, d, s,) we introduce the parameter 
Xi = rrii/T. For massless light quark-antiquark pairs one obtains &i(0) = 
2 • 9C(3)/2vr 2 w 1.10, 6 3 (0) = 2 • 7vr 2 /40 w 3.44 and 6 3 (0) = 1. For massive 
strange quark-antiquark pairs one can determine these factors exactly (see 
in Appendix). Here we investigate baryon symmetric systems, which means 
A u — A fl , Arf = Aj and X s = X s . Since boost invariance has been demonstrated 
to be a good approximation for the initially produced partons 0, we can 
estimate the initial parton fugacities, Ajj, A° and temperature To from 

1 dN 4 

^0 = ^2 T~ > £o = n -{k T ), (6) 

irR%n so dy n 
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where (kx) is the average transverse momentum. From HIJING one can obtain 
the initial values of the different gluon and quark fugacities. For different 
flavours the initial ratio is u : d : s ~ 2 : 2 : 1, approximately. 

Table 1. shows these relevant quantities at the moment r; so , for Au + Au 
collisions at RHIC and LHC energies. One can observe that the initial parton 
gas is rather hot reflecting the large average transverse momentum. However, 
the parton gas is undersaturated as compared to the ideal gas at the same 
temperature. The gas is also dominated by gluons with its quark content far 
below the chemical equilibrium value. We should emphasize that the initial 
conditions listed here result from HIJING calculation of parton production 
through semihard scatterings. Soft partons, e.g., due to parton production 
from the color field 



14 , are not included. 



Au+Au 


RHIC 


LHC 


T iso (im/c) 


0.7 


0.5 


e (GeV/fm 3 ) 


3.2 


40 


no(fm~ 3 ) 


2.15 


18 


(k±) (GeV) 


1.17 


1.76 


T (GeV) 


0.55 


0.82 




0.05 


0.124 




0.004 


0.01 


AS 


0.004 


0.01 


A° 


0.002 


0.005 



Table 1: Values of the relevant parameters characterizing the parton plasma 
at the moment r, so , when local isotropy of the momentum distribution is first 
reached. 



Master rate equations 

In general, chemical reactions among partons can be quite complicated because 
of the possibility of initial and final-state gluon radiations. Interference effects 
due to multiple scattering inside a dense medium, i.e., LPM suppression of 
soft gluon radiation has to be taken into account. One lesson one has learned 
from LPM effect [[L(J is that the radiation between two successive scatterings 
is the sum, on the amplitude level, of both the initial state radiation from the 
first scattering and the final state radiation from the second one. Since the 
off-shell parton is space-like in the first amplitude and time-like in the second, 
the picture of a time-like parton propagating inside a medium in the parton 
cascades simulations || shall break down. Instead, we shall here consider both 
initial and final state radiations together associated with a single scattering 
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(to the leading order, a single additional gluon is radiated, such as gg — * ggg), 
in which we can include LPM effect by a radiation suppression factor. The 
analysis of QCD LPM effect in Ref. 0] has been done for a fast parton 



traveling inside a parton plasma. We will use the results for radiation off 
thermal partons who average energy is about T, since we expect the same 
physics to happen. 

In order to permit the approach to chemical equilibrium, the reverse pro- 
cess, i.e., gluon absorption, has to be included as well, which is easily achieved 
making use of detailed balance. We consider only the dominant process 
gg — > ggg. Radiative processes involving quarks have substantially smaller 
cross sections in pQCD, and quarks are considerably less abundant than glu- 
ons in the initial phase of the chemical evolution of the parton gas. Here we are 
interested in understanding the basic mechanisms underlying the formation of 
a chemically equilibrated quark-gluon plasma, and the essential time-scales. 
We hence restrict our considerations to the dominant reaction mechanisms for 



the equilibration of each parton flavor. These are just four basic processes [15|: 



99 <-> 999, 99 <-> QiQi- (7) 

Other scattering processes ensure the maintenance of thermal equilibrium 
(99 99, 9Qi 9Q.ii e tc) or yield corrections to the dominant reaction 
rates (gqi <-> q { gg, etc.). 

Restricting to the reactions Eq. ([^D and assuming that elastic parton scat- 
tering is sufficiently rapid to maintain local thermal equilibrium, the evolution 
of the parton densities is governed by the master equations [^]: 

W) = (8) 

d^y) = 9,(%/) = ^ - ( 9 ) 

z V ' L qi' L g/ 

where h g = n g /X g and h qi = n q JXi denote the densities with unit fugacities, 
X g = X qi = 1, 03 and a l 2 are thermally averaged, velocity weighted cross 
sections, 

^3 = (<r(gg -> ggg)v), a\ = (a(gg -> qiqjv). (10) 

In Eq. (H) we neglected the exchange between quark flavours, because generally 
the gluonic channel is dominant. We have also assumed detailed balance and 
a baryon symmetric matter, n qi = n qi . If we neglect effects of viscosity due to 
elastic scattering |TJ, ITH] , we then have the hydrodynamic equation 



^(eO + P<9X = 0, (11) 
which determines the evolution of the energy density. 
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For a time scale r <C Ra, we can neglect the transverse expansion and 
consider the expansion of the parton plasma purely longitudinal, which leads 
to the Bjorken's scaling solution JlTj of the hydrodynamic equation: 

0. (12) 



de e + P 
dr r 



Evaluating this ultrarelativistic equation of motion by means of Eqs. 
one obtains 



X 9 J 14 
— + 4 1 

X g T t 3 



At T „ , . 1 
■^ + 4-D e {xi) + - 

Xi 1 T 



1 + g63(»i] 



^ = (13) 



where = e 9i + £g. and D e (0) = 63(0) = 1 for massless quark-antiquark pairs, 
D £ (x s ), b^(x s ) for massive strange quark can be found in Appendix. Rewriting 
the rate equations Eqs. (0)-@ i n terms of time evolution of the parameters 
T(r), X g (r) and Aj(r), one obtains 



A„ T 1 
X g T t 

K T 1 
^ + 3-D n ( Xl ) + - 

Xi 1 T 



R 3 (l-X g )-2j2RiU 



A? 



n n 



AL 



where the density weighted reaction rates R3 and i?2 are defined as 



Ra 



R% 



2 a 2 n g- 



(14) 
(15) 

(16) 



Notice that for a fully equilibrated system (X g = Aj = 1), Eq. (0) corresponds 
to the Bjorken solution, T(r)/T = (Vr) 1 / 3 . In Eq. (^) D n (0) = 1 and for 
D n (x s ) see Appendix. 



Parton equilibration rates 

To take into account of the LPM effect in the calculation of the reaction rate 
i?3 for gg — > (7 (7 (7 , we simply impose the LPM suppression of the gluon radia- 
tion whose effective formation time tqcd is much longer than the mean-free- 
path Xf of multiple scatterings to each gg — > ggg process. In the mean time, 
the LPM effect also regularizes the infrared divergency associated with QCD 
radiation. However, a 3 still contains infrared singularities in the gluon propa- 
gators. For an equilibrium system one can in principle apply the resummation 
technique developed by Braaten and Pisarski |18j to regularize the electric 



part of the propagators, though the magnetic sector still has to be determined 
by an unknown magnetic screening mass which can only be calculated non- 
perturbatively [TP| up to now. Since we are dealing with a non-equilibrium 
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system, Braaten and Pisarski's resummation may not be well denned. As an 



approximation, we will use the Debye screening mass [13 

& = % kf(k)dk = Ana s T 2 X g , (17) 
n z Jo 

to regularize all singularities in both the scattering cross sections and the 
radiation amplitude. 

The mean-free-path for elastic scatterings is 

(18) 



which depends very weekly on the gluon fugacity \ g [II 



The modified differential cross section for gg — > ggg has the form: 

Mlfe^*^i < ' (A '- TOOT)(,(Vi ^ lC08h! ' ) ' (19) 

where the second step-function accounts for energy conservation, s = 18T 2 is 
the average squared center-of-mass energy of two gluons in the thermal gas and 
the LPM suppression factor was approximated by a ^-function, 0(Xf — tqcd), 
where tqcd = cosh y/kj_ is the effective formation time of the gluon radiation 
in SU(3) flOf . We can complete part of these integrations: 

32 

Rs/T = —a s \ g (l + 87r« s A 9 /9) 2 X(A 9 ), (20) 

OCi\ 

where X(A 9 ) is a function of \ g , 

X{\ g ) = dx dz > v 



\ l + z ) 2 {x^j[x + (l + z)x D } 2 
1 cosh~ 1 (y / a;) 



sX f J [I + x(l + z)y D } 2 - Ax z y D 




+ 



where xd = fJ, 2 D X 2 , Vd — l^bl s an d the integral is evaluated numerically [TT|j . 

Note that in principle one should multiply the phase-space integral by 1/3! 
to take into account of the symmetrization of identical particles in gg — > ggg 
as in Ref. [BJ. However, for the dominant soft radiation we consider here, 
the radiated soft gluon does not overlap with the two incident gluons in the 
phase-space. Thus we only multiply the cross section by 1/2!. 

The quark equilibration rates R\ for gg — > can be also calculated, 
following the approximations of Ref. 0. We introduce the thermal quark 
mass as a cutoff in the logaritmically divergent integrals: 



M? = ml, + 



47T o . . 

—a s T 2 . (22) 
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Here m 0:i notes the bare quark mass, m 0iU = m ,d ~ 0, m , s = 150 MeV. 

The total cross section is dominated by the Compton diagrams. The cross 
section is the following for the process gg — > qiqf 



dai 
~dt 



3s 



log 



l + X 

i-x 



, IMf Mf\ ( 7 
1 + ^— --f )+X -7 



31 M 2 



2 s 



4 4 



(23) 



where x = y 1 — 4M 2 /s. 

To obtain the total cross section approximately, we integrate over thermal 
gluon distributions and we insert the thermal gluon mass, Mj, and the average 
thermal (s) = 18T 2 into the cross section: 



{(TiV) » - 



/C(2) N 


\ 2 rfcTj 




U(3), 




Mi] (s}=18T 2 



(24) 



By means of Eq. (24) one can derive numerically the production rates 
R\ = 1/2 ■ cr\n g which will be used in Eq. (|1^ -(|15|) in the calculation of the 
time evolution. 



Evolution of the parton plasma 

With the parton equilibration rates which in turn depend on the parton fu- 
gacity, we can solve the master equations self-consistently and obtain the time 
evolution of the temperature and the fugacities. The time dependence of T, X g 
and Aj are shown in Fig. l.a, and Fig. 2. a, for initial conditions listed in Table 
I. at RHIC and LHC energies. We find that the parton gas cools considerably 
faster than predicted by Bjorken's scaling solution ((T*) 3 r = const.) shown 
as dotted lines, because the production of additional partons approaching the 
chemical equilibrium state consumes an appreciable amount of energy. The 
accelerated cooling, in turn, slows down the chemical equilibration process, 
which is more apparent at RHIC than at LHC energies. Therefore, the parton 
system can hardly reach its equilibrium state before the effective temperature 
drops below T c w 200 MeV in a short period of time of 2-3 fm/c at RHIC en- 
ergy. At LHC energy, however, the parton gas will be closer to its equilibrium 
and the plasma may exist in a deconfmed phase for as long as 4-5 fm/c. 

We note that the initial conditions used here are the results of the HIJING 
model calculation in which only initial direct parton scatterings are taken into 
account. Due to the fact that HIJING is a QCD motivated phenomenological 
model, there are some uncertainties related to the initial parton production, as 
listed in Ref. 0. We can estimate the effect of the uncertainties in the initial 
conditions on the parton gas evolution by multiplying the initial energy and 
parton number densities at RHIC and LHC energies by a factor of 4. This 
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will increase the initial fugacities approximately by a factor of 4 at the same 
initial temperature. With these high initial densities, the parton plasma can 
evolve closer to an equilibrated gluon gas as shown in Fig. l.b, and Fig. 2.b, 
however the system is still non-equilibrated and dominated by gluons. The 
light quarks and strange quarks remain far from the full chemical equilibrium. 

Thus we can conclude that perturbative parton production and 
scatterings are very likely to produce a fully equilibrated parton 
plasma in ultrarelativistic heavy ion collisions at RHIC and LHC 
energies. The strangeness content will be far from the equilibrium, 
also. 



Gluon fragmentation 

We have seen that during the parton evolution strangeness and the other 
quarks were not produced at equilibrium level. Thus the last possibility to get 
some more strangeness is the gluon fragmentation before hadronization. In an 
oversimplified model we can investigate the relative strangeness enhancement 
or suppression in the hadronization. We assume a two step hadronization 
process at a certain critical temperature: (a) gluons decay into quark- ant iquark 
pairs forming quark-matter from the parton gas; (b) the quarks and antiquarks 
will be redistributed forming colorless hadrons. 

For simplicity, we will use the well-known string breaking factors to esti- 
mate the different flavour yields in the gluon decay process (g — > qi + gj, 
namely f u = fd = 0.425 and f s = 0.15. Generally, in a hot and under- 
saturated parton gas these factors will not be valid automatically, but in lack 
of correct, temperature and parton density dependent values, here we will 
use the above numbers. Let us define the relative strangeness content by the 
strange/non-strange flavour ratio g s : 

9s = (25) 
n u + n d 

One can determine easily this ratio g s in the non-equilibrated parton gas 
before and after fragmentation and we can compare its value to a fully equili- 
brated quark-gluon plasma characterized by the same temperature. Fig.l.c, d, 
and Fig.2.c, d, show the results of our calculation at RHIC and LHC energies. 

On Fig.l.c, full lines display the ratio g s for the previously calculated gluon 
dominated non-equilibrium parton gas (smaller value, g s ~ 0.25 and marked 
as "Non-eq.") and for the fully equilibrated quark-gluon plasma (larger value, 
g s 0.45 and marked as "Eq. "). Gluon fragmentation, characterized by above 
values fi, leads to the enhancement of the light quark flavours and the suppres- 
sion of the strangeness. On Fig.l.c, after fragmentation g s will drop in both 
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non-equilibrium (dotted line marked as "Fr-Non-eq" and g s ~ 0.18) and equi- 
librium cases (dashed line marked as "Fr-Eq" and g s ~ 0.28). Thus strangeness 
will be suppressed during hadronization in any case and any time step. Even 
if we consider 4 times higher initial parton fugacities (see Fig.l.d,) then we 
will have very much the same result as before. Moreover, the final g s ratio will 
depend very weakly on the initial condition. 

On Fig.2.c, we repeated the gluon fragmentation at LHC energy. We ob- 
tained an m 25% increase in strangeness comparing to the RHIC energy and 
an evolution which is slightly closer to the equilibrium ones. On Fig.2.d, the 4 
times larger initial parton fugacities display a similar increase in strangeness, 
however the strangeness ratio will also drop during hadronization. 

Our conclusions are the following: (a) before fragmentation g s is far from 
its equilibrium value and it depends very strongly on the initial parton fugac- 
ities; (b) the above simple gluon fragmentation will decrease this g s factor, 
suppressing the strangeness in a quark matter; (c) after fragmentation in equi- 
librium case g s ~ 0.3, in non-equilibrium case g s ~ 0.2 at RHIC energy and 
g s « 0.25 at LHC energy. 

This factor g s is important, because after gluon fragmentation the above 
flavour ratio will not change during the quark redistribution and it will deter- 
mine the final strange/non-strange hadron ratio. The hadron formation via 
quark and antiquark redistribution is beyond of the scope of our recent investi- 
gation. One can find a detailed description of such a hadronization mechanism 
in Ref. [|12|] where at 200 GeV/n bombarding energy (CERN SPS) it was ex- 
tracted g s = 0.16 for nucleon-nucleon and g s = 0.255 for S+S collision. 



Conclusions 

In this paper, we have calculated strangeness production in an equilibrating 
parton plasma, taking into account the evolution of the effective temperature 
and parton fugacities according to the solution of a set of rate equations. 
In the evaluation of the interaction rate for induced gluon radiation, a 
color dependent effective formation time was used which reduces the gluon 
equilibration rate through LPM suppression of soft gluons. 

We found that the thermal contribution during the parton equilibration 
with the current estimate of the initial parton density from HIJING Monte 
Carlo simulation are larger than the initial direct strangeness production, but 
the reached saturation level remains far from the full equilibrium at RHIC 
and LHC energies. The strangeness production depends on the initial condi- 
tion of the parton evolution. If uncertainties in the initial parton production 
can increase the initial parton density, e.g., by a factor of 4, the secondary 
strangeness production will be larger, but even in this case it will not reach 
the total chemical equilibrium. 
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However, during hadronization the yield of light quarks can be much larger 
as the strange one. Thus the hadronization process may lead to a strangeness 
suppression. This conclusion very much depends on the fragmentation ratios. 
The strange/non-strange quark ratio will be also far from its equilibrium value 
during the partonic evolution. Moreover it could drop during hadronization 
and the non-equilibrium characteristic will be conserved. 

Appendix 

For massive partons we have used the following expressions: 



bi(xi) 



, dj 
'2tt 2 

di 



2tt 2 



OO i 

3 ( 1 \n+l 1 



*?E(-i) 

n=l 

oo 



nxi 



o 1 4 
2— r • Xi 



£(-i) 

n=l 



n+l 



K 2 {nxi) 



3 1 

K 2 (nx i ) + - / -Ki(nxi) 



(nxi) 



(nxi)' 



i i 



3 (nxi) 



KAnXi 



(26) 
(27) 
(28) 



D n \Xi 



DJxi) 



E^i(-i) 



n+l 



KK 2 (nXi) + \Ki(nXi) 



(nxi) 



^{-lY^^Mnx,) 



E„° =1 (-i) n+1 


[ {n l )2 K 2 (nXi) + s ^K^nxi) + } 2 K (nXi)] 


E^=i(-l) n+1 


_ (n 1 Xi) 2K 2 (nx i ) + l (n 1 Xi) K 1 (nx i )_ 





(29) 
(30) 



with Xi = rrii/T and di = 6. The factor (— 1)™ +1 is connected to fermions. 
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